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ABSTRACT 

We suggest that the luminous extreme ultraviolet (EUV) emission which has been detected recently from clus- 
ters of galaxies is Inverse Compton (IC) scattering of Cosmic Microwave Background (CMB) radiation by low 
energy cosmic ray electrons in the intracluster medium. The cosmic ray electrons would have Lorentz factors of 
7 ~ 300, and would lose energy primarily by emitting EUV radiation. These particles have lifetimes comparable 
to the Hubble time; thus, the electrons might represent a relic population of cosmic rays produced by nonthermal 
activity over the history of the cluster. The IC model naturally explains the observed increase in the ratio of EUV 
to X-ray emission with radius in clusters. The required energy in cosmic ray electrons is typically 1-10% of the 
thermal energy content of the intracluster gas. We suggest that the cosmic ray electrons might have been produced 
by supernovae in galaxies, by radio galaxies, or by particle acceleration in intracluster shocks. 

Subject headings: cosmic rays — galaxies: clusters: general — intergalactic medium — radiation mechanisms: 
nonthermal — ultraviolet: general — X-rays: general 



1. INTRODUCTION 

Recently, observations of clusters of galaxies with the Ex- 
treme Ultraviolet Explorer (EUVE) Deep Survey instrument us- 
ing the Lexan B filter have revealed extreme ultraviolet (EUV) 
emission in excess of that expected from the previously ob- 
served. X-ray emitting intracluster medium (ICM). The pass- 
band of this instrument is 65-248 eV The excess EUV emission 
is seen in the Virgo cluster (Lieu et al. 1996b,c), Coma clus- 
ter (Lieu et al. 1996a), Abell 1795 (Mittaz, Lieu, & Lockman 
1997), Abell 2199, and Abell 4038 (Bowyer, Lieu, & Mittaz 
1997). The observed EUV luminosities in this band, corrected 
for Galactic absorption, range from Leuv 



; 9 X 10'*2/!-2 erg s" 



in Virgo to Leuv ~2x 10"*^/!^,, erg s"' in A1795 (where the 
Hubble constant is Ho = 5Qhso km s"' Mpc"'). In all of the 
cases except A2199, there is an associated soft X-ray excess 
seen in the ROSAT PSPC spectra of the clusters. The ratio of 
EUVE to PSPC fluxes and the correction for Galactic absorp- 
tion indicate that the emitted spectrum drops rapidly with pho- 
ton energy across the EUVE band. In A1795, the best-fit model 
for the spectrum implies that the average emitted photon energy 
is {hvEuv) ~ 75 eV 

The clusters with detected excess EUV emission are varied 
in their other properties. Virgo is not as rich as the others. Sev- 
eral have cooling flows in their centers, but Coma and A4038 
do not. Coma has a strong radio halo, but A 1795 and A2199 do 
not. Virgo is the closest cluster, while A 1795 is at a redshift of 
z = 0.0631. 

The EUV emission is extended over a scale comparable to the 
size of the cluster. In the three more distant clusters (A1795, 
A2199, and A4038), the EUV emission faUs off with radius 
more slowly than the cluster X-ray emission, so that the ratio of 



EUV emission to X-ray emission increases rapidly with radius. 
Virgo and Coma do not show such a trend, but are closer than 
the others, so the same angular scale doesn't correspond to as 
large a physical radius. 

In the papers reporting the discovery of this EUV emis- 
sion from clusters, a thermal emission model from diffuse gas 
was suggested. The temperature of this gas would be Tiian,, ^ 
(1 - 10) X 10^ K. While it is natural to suggest such a model, 
given that the X-ray emission from clusters is produced by ther- 
mal emission from diffuse gas with a higher temperature, there 
are a number of concerns with this hypothesis (Fabian 1996). 
First, aratherlargemassof this warm gas is required. InAI795, 
this mass is sa lO'^/jj^^M©, which is about the virial mass of 
the cluster (Mittaz et al. 1997). Moreover, gas at these temper- 
atures cools very rapidly. In A1795, the total cooling rate of 
gas would exceed 10^/zJqMq yr"'. There is no obvious energy 
source for this gas; Mittaz et al. suggest accretion of interclus- 
ter gas at a very high rate (^ IQr' Ii^-^Mq yr"'). In this thermal 
model, the bulk of the radiation would be line emission. A 
search for the O VI line, which occurs at longer wavelengths 
but is produced in gas of the same temperature, did not detect 
the line in some of the same clusters which show EUV emis- 
sion (Dixon, Hurwitz, & Ferguson 1996). Finally, the observed 
spatial distribution of EUV emission and the fact that the ratio 
of EUV to X-ray emission increases strongly with radius is not 
an obvious consequence of this model. 

Although a small amount of EUV emission might have been 
expected in the centers of cooling flow clusters if the gas contin- 
ues to cool below X-ray emitting temperatures, this cannot be 
the explanation of the observed EUV emission (Fabian 1996). 
First, the EUV emission is much too luminous, and would im- 
ply cooling rates ^10^ times higher than those seen in X-rays. 
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Second, not all of the clusters with EUV emission have cooling 
flows. Finally, the spatial distribution of the EUV emission is 
the opposite of that expected for cooling flow gas; in A 1795, for 
example, most of the EUV luminosity is at large radii, while the 
cooling flow gas is located at the cluster center where the gas 
density is highest. 

Fabian (1997) proposed a thermal model which avoids many 
of these problems. He suggested that the warm gas was located 
in turbulent mixing layers at the interface between embedded 
cold clouds (with T - lO'* K) and the ICM {T - 10** K). The 
power source for the EUV emission is then the thermal energy 
content of the ICM, which is large. Energetically, this model 
works for most of the clusters, but may have difficulties with the 
very luminous EUV emission of A 1795. The gas in the mixing 
layers is recycled, so that the cooling rates can be greatly re- 
duced. It is not clear that this model would predict the extended 
spatial distribution of the EUV emission when compared to the 
X-ray emission. 

2. INVERSE COMPTON EMISSION 

The difficulties with the thermal model for the EUV emis- 
sion have led us to consider a nonthermal model, in which this 
radiation is produced by Inverse Compton (IC) scattering of 
Cosmic Microwave Background (CMB) radiation by low en- 
ergy cosmic ray electrons. Because it has the largest luminosity 
and posed the greatest problems for the thermal model, we will 
use parameters for the EUV emission drawn from the obser- 
vations of A1795 (Mittaz et al. 1997). Consider a CMB pho- 
ton with a frequency vcmb which is scattered by a relativistic 
electron with an energy of -fnigC^, where 7 is the Lorentz fac- 
tor. On average, the resulting IC photon will have an energy 
of (hvic) = (4/3)7^ (/wcmb)- The typical photon in the CMB 
has an energy of {hvcMs) ~ ^.^kTcMB, where Tcmb = 2.73 K is 
the temperature of the CMB. If we associate hi'jc with the av- 
erage energy of EUV photons detected from a cluster, then the 
average Lorentz factor of the cosmic ray electrons is 



(7) w 300 



75 eV 



V2 



(1) 



Thus, the required cosmic ray electrons have very low energies 
compared to the typical energies of electrons which produce 
radio emission in radio sources. 

There is considerable evidence that the intracluster medium 
(ICM) also contains a magnetic field with a typical strength of 
B ^IfiG (e.g., Rephaeli, Gruber, & Rothschild 1987; Kim et 
al. 1990). The same low energy cosmic ray electrons will also 
generate synchrotron radio emission with an average photon 
frequency of (urati) = (55/96)(Vy7r)7^(eB/»vc)sina, where a 
is the pitch angle of the electron. This gives 



(l^rad) ;$ 5 X 10^ 



300 
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where the inequality comes from the fact that sina < 1. Thus, 
this radio emission would be well-below the Earth's iono- 
spheric cut-off, and unobservable. The ratio of the luminosity 
of the synchrotron emission to that of the IC emission is 
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where Ub and Uqmb are the energy density of the magnetic field 
and Cosmic Microwave Background, respectively. 

This shows that synchrotron radio emission produced by the 
same electrons would be unobservable. In particular, one would 
not expect that IC produced EUV emission would be restricted 
to clusters with radio halos. Such radio halos are, in fact, rather 
rare (e.g., Hanisch 1982). For example, A1795 has very strong 
EUV emission, but lacks a radio halo. Models in which the 
EUV emission from clusters is due to the electrons which pro- 
duce the radio halos have been proposed by Ensslin & Bier- 
mann (1997) and Hwang (1997). 

In addition to losses by synchrotron emission, it is important 
to consider if there are other energy loss processes for 7 « 300 
electrons in the ICM which would compete with IC scatter- 
ing of the CMB. First, the ICM contains other radiation fields 
which might also undergo IC scattering. These include the 
near-IR — optical light from galaxies, and the X-ray emission 
produced by the intracluster gas. However, the energy densities 
in these radiation fields are smaller than that in the CMB, and 
thus the losses are also smaller. Low energy cosmic ray elec- 
trons are also subject to losses from Coulomb interactions with 
thermal electrons in the ICM plasma and from bremsstrahlung. 
Both of these processes are proportional to the density of the 
ICM plasma. For electrons with 7 w 300 in an ionized hydro- 
gen plasma, the Coulomb losses are larger than those due to 
bremsstrahlung. The time scale for Coulomb energy losses is 
approximately (Rephaeli 1979) 
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V300/ 



10- 
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When compared with the time scale for energy loss by IC emis- 
sion [eq. (g) below], we find that IC emission dominates for 
plasma densities below a critical density ricnt. 
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10 



-3 (SL\ 
V300/ 



cm 



(5) 



This inequality is satisfied in most of the volume of a typi- 
cal cluster, and is particularly true at the outer radii where the 
EUV excess seems to be strongest (Mittaz et al. 1997). How- 
ever, at the centers of rich clusters, particularly in the cooling 
flows, the ICM density exceeds «„,> and Coulomb losses may 
be dominant. This might help to explain the lack of detectable 
EUV emission in these regions (Mittaz et al. 1997), although 
the brightnesses of the soft X-ray emission there also makes it 
more difficult to separate EUV and X-ray emission. 

IC losses will dominate for the low energy cosmic rays 
needed to produce the observed EUV emission from clusters. 
The lifetime of these particles is then 



tic=- 



"fMeC 



larcj^UcMB 



= 7.1 X 10 



V3OO/ 



yr. 



(6) 



Thus, the lifetime of a 7 « 300 electron is approximately the 
Hubble time, or the likely average age of clusters. We would 
then expect that any low energy cosmic rays with 7 « 300 
which were generated within or had escaped to the ICM over 
the lifetime of a cluster would still be present there. 

Note that equations (IT]) and (0) depend only on the properties 
of the CMB. Thus, these arguments show that the EUV spec- 
tral band is a unique region in which one can detect nonthermal 
emission from an accumulated relic population of cosmic rays 
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due to the total history of particle acceleration and activity over 
most of the lifetime of the universe in any diffuse regions (with 
fie <y ricrit)- To our knowledge, this unique role of the EUV 
spectral band has not been appreciated previously. On the other 
hand, the energy density of the CMB increases with redshift 
as UcMB fx (1 +z)'^. The lifetime for IC losses becomes much 
shorter than the Hubble time for redshifts z<,2. 



3. REQUIRED ENERGY IN COSMIC RAYS 

Let ncR('y)dj be the number density of cosmic ray electrons 
with Lorentz factors between 7 and -f + dj. Define an average 
Lorentz factor as 



(7) 



_ Jncrilh^d-y 



Jncrilhdj ' 
Then, the IC luminosity of the cluster is 
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Lie = T wUcmbEcr , 
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(7) 



(8) 



where Ecu is the total energy of the cosmic ray electrons in the 
cluster. If the Inverse Compton luminosity is equated to the ob- 
served EUV emission, the required total energy in cosmic rays 
electrons is found to be 



EcR w 2.4 X 10' 



62 



10**^ ergs 



300 



erg. 



(9) 



It is useful to compare the energy in cosmic rays to that in 
the thermal energy Egas of the intracluster gas. This ratio is 



EcR 



■^gas 
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10"*^ ergs 
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where Mgas is the mass of the ICM, and T is its average tem- 
perature. The ratio of pressures is lower by a factor of two, 
since the cosmic rays are relativistic and the thermal gas is not. 
Thus, the production of the observed EUV emission in clusters 
requires a relic cosmic ray electron population with an energy 
density and pressure which are 1-10% of that of the thermal 
gas. 

One concern is that the inclusion of cosmic ray ions may in- 
crease the total CR energy significantly. If the cosmic rays have 
a very steep energy spectrum, and if the number of cosmic ray 
ions at 150 MeV (7 w 300 in electrons) is reduced because these 
particles are not relativistic, then the increase in pressure due to 
ions might only be a factor of a few. On the other hand, cosmic 
ray ions outnumber electrons by a factor of ~10^ at the top of 
the Earth's atmosphere (e.g., Webber 1983). 



4. SPATIAL DISTRIBUTION OF EUV EMISSION 

The Inverse Compton model naturally explains the spatial 
distribution of the EUV emission and the fact the EUV to X- 
ray ratio increases rapidly with radius. The X-ray emission 
from the ICM is due to electron-ion collisions, and its emis- 
sivity is proportional to the square of the gas density. The In- 
verse Compton emission is due to collisions between cosmic 



ray electrons and CMB photons; the distribution of the latter is 
uniform. Thus, the Inverse Compton emissivity is linearly pro- 
portional to density, and the resulting surface brightness will 
decline with radius less rapidly than that for the X-ray emis- 
sion. 

As a specific but probably oversimplified model, let us as- 
sume that the pressure of the cosmic ray electrons, Pcr, is pro- 
portional to the thermal gas pressure Pgas- Further assume that 
the energy spectrum of the cosmic rays is independent of ra- 
dius. For this calculation, we will assume that the thermal gas 
pressure in a cluster is represented by the "beta model" with 
Pgas(r) oc [l + (r/rc)]~^^^. Here, r is the radius and r^ is the core 
radius. Then, the EUV surface brightness varies as lEuvir) c)C 
[l+(r/rc)]~^^^'^^^, while the X-ray surface brightness varies as 
Ix(r) oc [1 +(r/rc)]-3^+'/2 'phe ratio of the EUV to X-ray surface 
brightness goes as Ieuv/Ix oc VPgas oc [l+ir/rc)]^^'^. 




radius (arcniin) 

Fig. 1. — The ratio of EUV to X-ray surface brightness in A1795 as a 
function of projected radius (Mittaz et al. 1997). The solid curve is the "beta 
model" prediction for this ratio in the Inverse Compton theory, assuming a con- 
stant cosmic ray to gas pressure ratio. The parameters of the beta model are 
taken from the fit to the X-ray surface brightness by Briel & Henry (1996). 

Figure [I] shows the observed ratio of EUV to X-ray surface 
brightness measured for several annuli in A 1795 (Mittaz et al. 
1997). For comparison, the results of this simple beta model 
are shown. The parameters of the beta model (jc and /3) were 
taken from the fit to the observed X-ray emission by Briel & 
Henry (1996); there are no free parameters in this model ex- 
cept for the overall normalization. Given the many simplifying 
assumptions in this model, the detailed agreement is probably 
fortuitous. However, the comparison does show that the Inverse 
Compton model can explain the qualitative radial behavior of 
the EUV surface brightness. 

5. DISCUSSION 

Normally, nonthermal emission processes produce broad 
power-law spectra, partly as a result of power-law energy spec- 
trum of the cosmic ray electrons. The observations of clus- 
ters suggest that the EUV emission drops rapidly in going from 
EUV to soft X-ray photon energies. On the low energy side, the 
total energies involved would be large and the surface bright- 
nesses would violate observational limits if the spectrum con- 
tinued to rise down to the near UV and optical bands. However, 
there are reasons to think that the CR electron population in 
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clusters might peak near the energies required to produce the 
EUV emission (7 « 300). On the higher energy side, the elec- 
tron population will be depleted by IC losses over the age of 
the cluster [eq. (^)]. For continual particle injection, this steep- 
ens the exponent of a power-law energy spectrum by unity. If 
most of the CRs were injected early in the history of the cluster, 
the effect would be larger. On the lower energy side of the CR 
spectrum, electrons will lose energy rapidly to Coulomb losses 
even at low ICM densities [eq. (Q)]. 

What is the ultimate source of these relic cosmic ray elec- 
trons, which have accumulated over the lifetime of the cluster? 
Note that we typically need Ecr '^ 3 x 10*"' erg, although the 
values in A 1795 are a factor of 10 larger [eq. (ph]. This im- 
plies that the ratio of the cosmic ray electron energy to that of 
the thermal gas is typically ^1%, but with a value of ~10% in 
A1795. 

One possible source for these cosmic rays would be super- 
novae, supernova remnants, and pulsars associated with the 
stellar population in the galaxies. The present rate of star for- 
mation in the early-type galaxies in rich clusters is rather low, 
and the Type la supernovae produced by the older stellar pop- 
ulation occur at too low a rate to give the required number of 
cosmic rays. However, it is believed that these galaxies had 
much higher supernova rates in the past. The heavy element 
abundances in clusters provide a useful constraint on the over- 
all number of such supernovae. As a simple example, let us as- 
sume that each supernova produces 0.5 Mq of heavy elements 
and a total cosmic ray electron energy (from the supernova, its 
remnant, pulsar, etc.) of 10**^ erg. Then, the production of a 
half-solar heavy element abundance in the ICM will generate a 
cosmic ray electron population of Eci/Egas ^ 1%. This might 
be barely adequate. One concern is that the cosmic ray elec- 
trons might experience strong Coulomb or adiabatic losses be- 
fore they reach the ICM. Another worry is that the early-type 
galaxies in clusters may have experienced most of their star 
formation and supernovae before z ^ 2, and the cosmic rays 
would have also undergone strong IC losses. A prediction of 
this model for the origin is that the EUV properties of clusters 
might be relatively uniform from cluster to cluster, just as the 
abundances of heavy elements in the ICM are fairly constant. 

Another possibility is that the cosmic rays may have come 
from radio galaxies, quasars, or other AGNs in the cluster over 
its lifetime. The total cosmic ray electron energy required cor- 
responds to 10-100 Cyg-A radio sources, and the present day 
radio luminosity function of clusters suggests that there are too 



few bright radio galaxies to produce the cosmic rays. However, 
the energy content of radio sources and their radio luminosities 
are not necessarily well-correlated. Moreover, there is evidence 
that both the general population of luminous radio galaxies and 
quasars was higher in the past, and that they were more com- 
monly associated with clusters of galaxies. In this model, one 
might expect large variations in the EUV luminosity of other- 
wise similar clusters, which would be associated with the small 
number statistics of particularly energetic AGN. 

Finally, the cosmic ray electrons might have been gener- 
ated in situ in the ICM by particle acceleration associated with 
shocks and/or strong turbulence. Unless the ICM underwent 
strong preheating, most of the gas must have passed through 
high Mach number shocks with velocities ~ 10^ km s"'. These 
would include overall cluster accretion shocks (if much of the 
cluster formed from a single large perturbation) and/or subclus- 
ter merger shocks if the cluster formed hierarchically. In the 
interstellar medium of our Galaxy, shocks of this velocity asso- 
ciated with supernova remnants also produce cosmic ray elec- 
trons with about 1-10% of the shock energy. This could lead to 
a population of such electrons in the ICM with Eci/Egas ~ 1- 
10%. In this model, modest variations in the EUV luminosity 
from cluster to cluster might be expected as a result of varia- 
tions in their detailed histories. 

The IC model for the EUV emission of clusters predicts a 
continuum spectrum for the emission. On the other hand, in 
the thermal model most of the emission is in lines. Thus, 
spectra of this component could decide between thermal and 
nonthermal models. The cosmic ray electrons required by the 
IC model would also generate some gamma-ray emission via 
the scattering of higher energy photons and bremsstrahlung. 
The expected fluxes are 10"'^- 10~" erg cm"^ s"' at 0.3 MeV 
(from IC scattering of optical photons) and at 150 MeV (from 
bremsstrahlung). These predicted fluxes are near the sensitivity 
limits with current instrumentation, but should be detectable in 
the future. 
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